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Abstract. We report a phonon density of states measurement of a-ZnCLo using the 
coherent inelastic neutron scattering technique and a lattice dynamical calculation 
in four crystalline phases of Z11CI2 using a transferable interatomic potential. The 
model calculations agree reasonably well with the available experimental data on the 
structures, specific heat, Raman frequencies and their pressure variation in various 
crystalline phases. The calculated results have been able to provide a fair description 
of the vibrational as well as the thermodynamic properties of ZnCl2 in all its four 
phases. 

PACS numbers: 63.20.Dj, 78.70.Nx, 64.30. +t, 65.40.De 
1. Introduction 

Among the metal halides of form MX 2 (M: metal, X: halide), Z11CI2 continues to draw 
attention in current research for several reasons. Firstly, the crystalline ZnCl2 has four 
structurally determined polymorphs^ |2] at the ambient pressure, viz. a-ZnCl 2 (I42d, 
Z=4], /3-ZnCl 2 {P2 x /n, Z=12), 7-ZnCl 2 (P4 2 /nmc, Z=2), and <5-ZnCl 2 (Pna2 1 , Z=4) in 
which the Zn atom is tetrahedrally coordinated with CI atoms. Secondly, ZnCl 2 displays 
unique phase transitions upon compression as well as pressure release. For instance, 
the four-coordinated a-ZnCl 2 transforms into a six-coordinated CdCl 2 -like phase at 
2.1 GPa, which further converts into four-coordinated 7-ZnCl 2 upon decompression [3]. 
Thirdly, molten ZnCl 2 , which is commercially used as a solvent and catalyst in the hydro 
cracking of coal slurries and heavy crude oil fractions jl], typifies high viscosity EJ EI) 
low ionic diffusivitiesjH], and very low electric conductivity [H] due mainly to its non- 
uniform polymeric network linkage. [TU] Further, liquid ZnCl 2 has the ability to supercool 
into a glass at a relatively low glass transition temperature (T g = 375 K). Finally, 
ZnCl 2 glass, which is intermediate between strong (e.g. Si0 2 ) and fragile (e.g. CKN) 
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glasses, finds its potential use as a low-loss optical transmission medium in the near 
infrared region. [T2*] 

All the crystalline phases of ZnCl2 at ambient pressure are formed by the network 
of ZnCU tetrahedra sharing their corners. But this network is rather soft (unlike SiC>2 or 
BeF 2 ) resulting in a strongly hybridized character for the normal modes of vibration. [TT] 
Further, in the low frequency range (< 20 cm -1 ), where the collective acoustic modes 
are expected to dominate the vibrational density of states, liquid ZnCl 2 exhibits a 
prominent boson peak|13[ HI] reflecting its more complex microstructure. Ribero et 
al.|15j examined the nature of instantaneous normal modes in the atomic motion of 
liquid ZnCl2, which involve the coupled vibration of several quasi- molecular ZnCl4 units 
of the network and has the local torsional character. Although molecular dynamics 

(MD) and Monte Carlo (MC) simulations have been reported for glassy as well as 
molten ZnCl2, the need for a realistic interatomic potential has remained unfulfilled. |16j 
The present work is primarily intended to make an attempt in formulating an effective 
transferable interatomic potential that can be applied to any of the polymorphs of ZnCl2 
with equal ease. Angell et al.jTT] had already suggested that there is little distinction 
between the glassy and the crystalline behavior of ZnCl2 so far as the dynamics is 
concerned. They further opined that the short time dynamics of liquid ZnCl2 near 
its melting point (T m = 593 K) might be regarded as solid-like in character. Later 
on, X-ray diffraction measurements of molten ZnCl2 by Triolo and Narten^H] favored 
this argument strongly. Anomalously low enthalpy and volume changes on fusion are 
supposed to be behind the resemblance of ZnCLj melt with the crystalline structure. [I] 
An extensive neutron diffraction study by Desa et al.^H] further showed that vitreous 
ZnCl2 has a distorted closed packed array of CI atoms in which Zn atoms are randomly 
distributed in the tetrahedral sites. Numerous such studies undoubtedly imply that if 
we are able to understand the dynamical properties of crystalline polymorphs of ZnCl2, 
microscopic interpretation for glassy and molten phases will also be robust. In the 
present work we have also determined the phonon density of states in the a-phase by 
coherent inelastic neutron scattering. The potential model calculation is found to be in 
fair agreement with the neutron data. 

2. Crystal structure 

Till 1961, it was commonly believed that the ambient stable phase of ZnCl2 had a layered 
structure in which each Zn atom was coordinated to six CI atoms to form regular 
octahedra with van der Waals bonding between the layers. [2*U] Later investigations 
proved that ZnCL> has bonding based on a tetrahedral framework. Initially three solid 
phases of ZnCl2 were known crystallographically. These were a, (3 and 7 polymorphs of 
ZnCl2 belonging to monoclinic and tetragonal crystal classes (Table H}. While a and /3 
phases have cristobalite type three dimensional networks of corner linked tetrahedra P, 
7-ZnCl2 crystallizes to a layered structure of cross linking tetrahedra as in case of 
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red HgI 2 .j2I] However, in 1978, Brynestad and Yakel[2*| 122] pointed out that all three 
modifications are rather the result of water contamination and that purely anhydrous 
ZnCl2 has an orthorhombic unit cell containing corner linked ZnCLi tetrahedra. This 
orthorhombic phase has been referred to as £-ZnCl2 in various papers. [3 EI] Fig. E 
suggests that the close packed layers in <5-ZnCi2 are horizontal. Further, temperature 
has a significant impact [22] on the structure of ZnCl2. While rapid cooling of the melt 
(or devitrification of the glass) produces the a-phase, slow cooling of the melt favors 
the formation of /3-ZnCl2. On the other hand, concentrated solution yields the 7-phase 
upon crystallization. 

3. Experiments and data analysis 

Experimental studies involving crystalline ZnCi2 encounter difficulties because of its 
extreme hygroscopicity. The reagent grade chemical ZnCl2 (purity 98%), purchased 
from Merck Chemical Co. (Germany), was dried for a period of 8 hours. It was then 
exposed into the glove box under dry argon atmosphere for transferring into different 
containers. About 30 g of the polycrystalline ZnCl2 sample so obtained was placed in a 
thin walled aluminum sample holder of circular shape for the inelastic neutron scattering 
(INS) studies. Another cylindrical vanadium can was made use of for taking the neutron 
diffraction pattern of the same sample. Neutron powder diffraction confirmed that 
the sample belongs to the a-phase. The structural refinements of the powder data 
were carried out using the Rietveld profile refinement method with the program DBWS 
9411. [21] The cell constants and atomic coordinates refined from the data recorded at 
300 K are in close agreement with the reported data, 

Inelastic neutron scattering measurements were carried out at 300 K using a medium 
resolution triple- axis spectrometer |25j at the Dhruva Reactor, Trombay. The instrument 
used a monochromated beam of neutrons produced by diffraction from the (111) planes 
of a copper crystal. Neutrons scattered inelastically by the sample were again diffracted 
by a pyrolytic graphite (002) analyzer and subsequently detected in a 10 BF 3 proportional 
counter. All the measurements were made in the energy loss mode with constant 
momentum transfer (Q). The elastic energy resolution was about 15% of the initial 
energy. Several scans were performed with the final energy (Ef) values of 30 meV and 
Q values of 5 to 6 A~ x that are much larger than the size of the Brillouin zone (about 1 
A -1 ). Additional measurements of the neutron background were made by detuning the 
analyzer crystal from the position by ±5° and repeating the scan for the same sample. 

The neutron weighted phonon density of states, g n (E), are obtained from the 
measured scattering function S(Q, E) through the following relation [2T>] 




(1) 



(2) 
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where n(E,T) = [exp(E/KT) — l] -1 ; A, B are normalization constants; 2W(Q) is 
the Debye- Waller factor; b p , M p , and g p {E) refer respectively to the neutron scattering 
length, mass and partial phonon density of states of the p th atom in the unit cell. The 
quantity within (...) represents the average over all Q values. The factor ^ £ turns out 
to be 0.063 and 0.474 barn/a. m.u. for Zn and CI atoms respectively. 



4. Model interatomic potential 



As previously stated, there is a need to formulate a realistic potential model that can 
effectively model the structure and dynamics in various polymorphs of ZnCl2. In the 
present model, Born-Mayer repulsive and attractive van der Waals interactions among 
different pairs of atoms have been incorporated. Because of the covalent nature of the 
Zn-Cl bond, a stretching term has to be included further in the potential model. The 
model thus turns out to be 



V(r) = A exp 



Br 



a 



D exp 



n{r — r^y 



2r 



(3) 



R(k) + R(k' 

where A (=1822 eV) and B (=12.364) are two constants, used extensively in several 
previous works. fI7 \ I28 | R(k) and R(k') refer to effective radius parameters for Zn 
and CI atoms. Ckk' accounts for van der Waals terms between Zn-Zn, Zn-Cl and Cl-Cl 
atoms; ro is the equilibrium Zn-Cl bond length; D and n are two adjustable parameters. 
At the observed lattice constants and atomic positions in the unit cell, the parameters 
were so adjusted as to give nearly a zero internal stress and net force on each individual 
atom. Once the structural constraints are met, the potential parameters were further 
adjusted so that the calculated eigen frequencies are real for all the wave vectors in 
the entire Brillouin zone. The group theoretical analysis provides the symmetry vectors 
necessary for block diagonalization of the dynamical matrix. Table El lists the final 
values of the optimized parameters. All the necessary calculations were made using 
DISPRjHn] developed at Trombay. Lattice parameters and atomic coordinates obtained 
as a result of the potential minimization for various polymorphs of ZnCl2 are displayed 
in Tables El IH El & El At high pressures, the equilibrium crystal structures were obtained 
from the free energy minimization. 



5. Vibrational modes: Raman and infrared active 



A group theoretical classification of various zone-center phonon modes for each of the 
crystalline modifications in ZnCl2 is given in Table ^ Experimental mode assignments 
are available for only two polymorphs of ZnCl2 (viz. a-ZnCl2 and 7-ZnCl2). Measured 
optical phonon data|21[IHll221EIj are found to be in good agreement with our calculated 
results (Tables El El & fTU|) . It may be noted that the lowest frequency mode in the 
7-phase (36 cm -1 ) is almost half the frequency of that in the a-phase (80 cm -1 ). Further 
the pressure dependence of Raman modes in a, 7 and 5 forms of ZnCl2 are calculated. 
Except for the a-phase, in all the others, the calculated pressure dependence compares 
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well with the available measured PI 13*2*] data (see Fig- EJ) - Since /3-ZnCl 2 has as many as 
108 phonon modes, pressure dependence of these modes is not shown here. 

6. Phonon dispersion relations and phonon density of states in ct-ZnCl 2 

As Table "*0 suggests, there is a total of 18 degrees of freedom in a-ZnCl 2 , giving rise 
to 18 phonon frequencies. Following are the representations of normal modes along the 
three high symmetry directions (viz. E, A and A). 

£(100): 8Ei + 10£ 2 
A(001): 4Ai + 4A 2 + 5A 3 (A 3 being doubly degenerate) 
A(110): 9Ai + 9A 2 

Fig. El displays the calculated phonon dispersion relations in a-ZnCl 2 . Considerable 
dispersion and anti-crossing are found for the branch in all the three symmetry directions 
which indicate strong hybridization. There are no indications of any phonon softening 
in the calculations in any of the four polymorphs at ambient pressure. 

We have plotted the measured phonon density of states along with the calculated 
spectrum in Fig. EJ The experimental one-phonon spectra are obtained by subtracting 
the multiphonon contribution from the experimental data. The multiphonon 
contribution is obtained using the Sjolander formalism. [*"3~] The experimental data 
(Fig. [""J confirm two broad structures with a band gap around 25 meV. The measured 
neutron data[3"""j f° r vitreous ZnCl 2 and the Raman density of states data|2H] for molten 
ZnCl 2 are qualitatively similar to the present neutron data on crystalline a-ZnCl 2 . The 

computed partial phonon density of states is shown in Fig 5(a). This indicates that the 
Zn ions contribute significantly at energies below 10 meV while from 10 meV to 19 meV, 
the vibrations of CI ions dominate. However, in the higher energy region, contributions 
from both the ions remain more or less similar. A comparison of the calculated phonon 
density of states for all the crystalline polymorphs is shown in Fig. 5(b). Two band 
gaps (a broad one around 23 meV and a very narrow one around 36 meV) are observed 
in each of the phases, though of slightly varying degree. For these calculations, all the 
phonon modes are integrated at each wave vector within the irreducible Brillouin zone 
on a 10 x 10 x 10 mesh with an energy resolution of 1.0 meV. 

7. Thermodynamic properties 

7.1. Heat capacity and Debye temperature 

As Fig. EKa) suggests, there is a fairly good agreement between the calculated and the 
experimental data]?""""] on the molar heat capacity. As the inset shows, there is vary small 
variation among heat capacities of various polymorphs, which reflects the small variation 
in the respective phonon density of states. In particular, the spike in AC P at about 50 
K in the inset of Fig. H""(a) arises due to the different phonon density of states of the 
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a-phase at low energy in Fig. \5[h). Fig. E(b) displays the Debye temperature {6d) plot 
for different crystalline modifications in ZnCl2. The characteristic Debye temperature 
in the present model turns out to be 375 K, while the experimental value is 340 K|36j. 

7.2. Mode Griineisen parameters and thermal expansion 

The mode Griineisen parameters are important indicators of the anharmonicity of 
phonon modes and their role in determining the thermal expansion. For the i th mode 
the Griineisen parameter is given by 

din uji , A . 

In the quasiharmonic approximation, the zth phonon mode contributes (^y)j%Cyi 
to the volume thermal expansion (ay), where Cvi denotes the contribution of the i th 
mode to the volume heat capacity. As heat capacity and bulk modulus values remain 
similar for each of the polymorphs, it is obviously the mode Griineisen parameter which 
dominates in the determination of ay. Fig. [7| displays averaged ji as a function of 
phonon energy for different crystalline polymorphs in ZnCl2. One may notice that while 
for a, (3 and S phases, 18 meV phonons have maximum variations with volume, 7-ZnCl2 
behaves slightly differently. This is further reflected in the subsequent thermal expansion 
calculation (see Fig. EJ). 



7. 3. Equation of state 

The crystal structure parameters as a function of pressure (at T=0 K) are calculated 
corresponding to the minimum of Gibbs free energy. Fig.[H]shows how the cell dimensions 
belonging to different crystalline polymorphs of ZnCl2 vary with pressure. It is observed 
that the compressibility is anisotropic. While c axis is the most compressible in (3 and 
7 forms of ZnCl2, it is the least compressible in a-ZnCl2. In 5-ZnCl2, b axis is found 
to be the most compressible. However, the volume derivative with pressure in all the 
phases looks a bit alike, as their bulk moduli are similar. 



8. Conclusions 



We have developed an interatomic potential model for ZnCl2 that mimicks the 
vibrational and thermodynamic properties of the polymorphs of crystalline ZnCl2 and 
reproduces fairly well the inelastic neutron scattering data for the a-phase. The fact 
that the 7-phase shows differences in Griineisen parameters and lattice expansion is 
essentially due to its layered structure, which is very different from other corner-shared 
polymorphs. This transferable potential model may be exploited further to investigate 
the properties of the glassy and the liquid phases of ZnCl2. 
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Figure 1. Crystal structures of polymorphs of ZnCl2, viz. (a) a-ZnCh, (b) /3-ZnCi2, 

(c) 7-ZnCl2 and (d) i5-ZnCl2. The small circles represent CI atoms, while Zn atoms 
are located inside (not visible here) the ZnCU tetrahedra. 

Figure 2. Pressure dependence of phonon modes in (a) a-ZnCl2 (b) 7-ZnCl2 and (c) 
<S-ZnCl 2 . 

Figure 3. Calculated phonon dispersion relations in a-ZnCl2. 

Figure 4. Plots of (a) experimental and (b) calculated neutron weighted phonon 
density of states for a-ZnCl2. The multiphonon contribution at 300 K has been 
subtracted from the experimental data to obtain the experimental one-phonon 
spectrum. The corresponding calculated spectrum has been convoluted with the energy 
resolution of the instrument. 

Figure 5. (a) Calculated partial phonon density of states of various atoms in a-ZnCl2. 
(b) A comparison of calculated one phonon density of states in various crystalline 
polymorphs of ZnCi2. 

Figure 6. Plots of (a) heat capacity (Cp) and (b) Debye temperature {9d) for all the 
crystalline polymorphs of Z11CI2. The inset in (a) shows the specific heat differences 
of various polymorphs with respect to the a-phase. 

Figure 7. Calculated plots of the mode Griineisen parameter and the volume thermal 
expansion coefficients of ZnCl2 in its various crystalline phases. 

Figure 8. Calculated equations of state for crystalline ZnCl2 in (a) a (b) (3 (c) 7 and 

(d) 8 phases. 



Collective dynamics in crystalline polymorphs of ZuCIq. 
Tables and table captions 



10 



Table 1. Group theoretical classification of zone-centre phonon modes in various 
crystalline polymorphs of ZnCl2. Raman and infra-red active phonon branches are 
indicated within the bracket by R and ir respectively. 



Phase Crystal system Space-group Z Group theoretical classification at zone-centre 



a 


Tetragonal 


I42d 


4 


A x (R) + 2A 2 + 2B X (R) + 3B 2 (R, ir) + 5E(R, ir) 


P 


Monoclinic 


P2 1 /n 


12 


27 Ag(R) + 27 A u {ir) + 27 B g {R) + 27 B u (ir) 


7 


Tetragonal 


P4: 2 /nmc 


2 


Aig(R) + 2A 2u (ir) + 2B lg (R) + B 2u (ir) + 3E g (R) + 3E u (ir 


5 


Orthorhombic 


Pna2 1 


4 


9Ai(R, ir) + 9A 2 (R) + §B X {R, ir) + 9B 2 (R, ir) 



Table 2. Transferable potential parameters for ZnCl2 polymorphs. 



R(k) 


Ci^eVA 6 ) 


D{eV) 


nil- 1 ) r (A) 


Zn CI 


Zn-Zn Zn-Cl Cl-Cl 






1.03 2.52 


20 135 550 


3.00 


4.60 2.30 



Table 3. Lattice constants and position parameters for a-ZnCl2. For the space group 
IA2d : Zn atoms are located at (0,0,0) and (^,0, |) while CI atoms are located at 
general positions of (x, \, |), (x, \, \), {\,x, |), |) respectively. 

This work Experimental data 
(Ref. 1) 

a(A) 5.406 5.410 
c(A) 10.415 10.300 
V(A 3 ) 304.34 301.46 
x 0.24 0.25 
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Table 4. Lattice constants and position parameters for /J-ZnCLj. For the 
space group P2i/n, Zn and CI atoms are located at general positions of (x,y,z), 
(x + + \,z + |), (x,y,z), (x + |, y + |, 2 + |) respectively. 



This work Experimental data 

(Ref. 1) 



afA) 




6.519 






6.500 




6fA) 




11.330 






11.300 














1 9 onn 




/?(") 




89.8 






90.0 




v(A 3 ) 




906.43 






903.43 




Atom 


X 


y 


z 


X 


y 


z 


Zn 


0.165 


0.167 


0.065 


0.167 


0.167 


0.063 


Zn 


0.166 


0.499 


0.186 


0.167 


0.500 


0.188 


Zn 


0.666 


0.667 


0.187 


0.667 


0.667 


0.188 


CI 


0.328 


0.003 


0.127 


0.333 


0.000 


0.125 


CI 


0.326 


0.334 


0.124 


0.333 


0.333 


0.125 


CI 


0.335 


0.663 


0.128 


0.333 


0.667 


0.125 


CI 


0.836 


0.167 


0.127 


0.833 


0.167 


0.125 


CI 


0.835 


0.503 


0.127 


0.833 


0.500 


0.125 


CI 


0.828 


0.831 


0.122 


0.833 


0.833 


0.125 



Table 5. Lattice constants and position parameters for 7-Z11CI2. For the space group 
Pi2/nmc, Zn atoms are located at (0,0,0) and (^, ^, ^) while CI atoms are located at 
general positions of (0, \,z), (0, \, z + i), (i, 0, z + |), (5, 0, zj respectively. 





This work 


Experimental data 






(Ref. 1) 


a(A) 


3.733 


3.700 


c(A) 


10.859 


10.670 


v(A 3 ) 


151.31 


146.07 


z 


0.121 


0.125 
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Table 6. Lattice constants and position parameters for <5-ZiiCl2- For the space group 
Pna2\ 1 Zn and CI atoms are located at general positions of (x,y,z), (x, y, z + h), 
(x + |, y + i, z) , (x + i, y + h, z + i) respectively. U(A 2 ) gives the isotropic thermal 
amplitudes. 



This work Experimental data 

(Ref. 2) 



a(A) 




6.444 






6.443 




6(A) 




7.666 






7.693 




c(A) 




6.122 






6.125 




/?(°) 




90.0 






90.0 




v(A 3 ) 




302.40 






303.59 




Atom 


,r 


y z 


U 


X 


y z 


U 


Zn 


0.081 


0.125 0.378 


0.0303 


0.0818 


0.1251 0.3750 


0.0253 


CI 


0.075 


0.122 0.005 


0.0296 


0.0702 


0.1223 0.0041 


0.0305 


CI 


0.084 


0.630 0.003 


0.0296 


0.0841 


0.6332 0.0062 


0.0292 



Table 7. Comparison of calculation and experimental data of optical modes in a- 
ZnCl2- Modes of A2 representation are not optically active. 

Optical phonon modes (cm -1 ) 
Representations This work Experimental data 

(Ref. 3) (Ref. 23) (Ref. 31) 





223 


226 


233 


245 


A 2 


147 
301 








Bx 


94 
273 


117 




113 


D, 


111 

325 


128 




128 


E 


75 


76 


82 


80 




99 


100 


103 


103 




263 










327 
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Table 8. Comparison of calculation and experimental data of optical modes in (3- 
ZnCl 2 . 

Optical phonon modes {cm" 1 ) 

This work Experimental data 

(Ref. 31) (Ref. 23) 



4 




4 


R 


(A B ) (A 


34 


33 


34 


36 




40 


47 


43 


52 


42 


44 


49 


48 


57 


57 


47 


55 


53 


61 


61 


48 


(il 


58 


69 


70 


57 


66 


G4 


74 




G4 


71 


70 
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Table 9. Comparison of calculation and experimental data of optical modes in 7- 
ZnCl 2 . 



Optical phonon modes (cm 1 ) 
Representations This work Experimental data 

(Ref. 3) (Ref. 31) 



Alg 


244 
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A 2u 


258 
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66 
271 




B 2 u 


216 




E g 


44 


36 38 
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E u 


68 
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Table 10. Comparison of calculation and experimental data of optical modes in 6- 
ZnCl 2 . 

Optical phonon modes (cm" 1 ) 





This work 




Experimental data 










(Ref. 3) 


(Ref. 16) 
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